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Abstract 

In this paper, we present the diversity order analysis of bit-interleaved coded multiple beamforming (BICMB) 
combined with the constellation precoding scheme. Multiple beamforming is realized by singular value decompo- 
sition of the channel matrix which is assumed to be perfectly known to the transmitter as well as the receiver. 
Previously, BICMB is known to have a diversity order bound related with the product of the code rate and the 
number of parallel subchannels, losing the full diversity order in some cases. In this paper, we show that BICMB 
combined with the constellation precoder and maximum hkeUhood detection achieves the full diversity order. We 
also provide simulation results that match the analysis. 

I. Introduction 

When the perfect channel state information is available at the transmitter to achieve spatial multiplexing 
and thereby increase the data rate, or to enhance the performance of a multi-input multi-output (MIMO) 
system, beamforming can be employed |[T1. The beamforming vectors are designed in [|2l, lO for various 
design criteria, and can be obtained by singular value decomposition (SVD), leading to a channel- 
diagonalizing structure optimum in minimizing the average bit error rate (BER) Q. 

It is known that an SVD subchannel with larger singular value provides larger diversity gain. During 
the simultaneous parallel transmission of the symbols on the diagonalized subchannels, the performance 
is dominated by the subchannel with the smallest singular value, resulting in losing the full diversity order 
BH, [|5l. To overcome the degradation of the diversity order of multiple beamforming, bit-interleaved coded 
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multiple beamforming (BICMB) was proposed |l6l, Q. This scheme interleaves the codewords through 
the multiple subchannels with different diversity orders, resulting in better diversity order. BICMB can 
achieve the full diversity order offered by the channel as long as the code rate Rc and the number of 
subchannels used S satisfy the condition R^S < 1 [HI. 

We showed in [9] and [10] that constellation precoded multiple beamforming, which converts a symbol 
into a precoded symbol and distributes the precoded symbol over the subchannels, can compensate for 
the diversity loss caused by the uncoded multiple beamforming. In this paper, by calculating pairwise 
error probability (PEP), we present the diversity analysis of Bit-Interleaved Coded Multiple Beamforming 
with Constellation Precoding (BICMB-CP), which adds the constellation precoding stage to BICMB. We 
show that adding the constellation precoder to the BICMB system which does not satisfy the full diversity 
condition guarantees the full diversity order when the subchannels to transmit the precoded symbols are 
properly chosen. Simulation results are shown to prove the analysis. 

The rest of this paper is organized as follows. The description of BICMB-CP is given in Section [III 
Section |in] presents the diversity analysis through the calculation of the upper bound to PEP. Simulation 
results supporting the analysis are shown in Section |IVl Finally, we end the paper with our conclusion in 
Section |Vl 

Notation: Bold lower (upper) case letters denote vectors (matrices). diag[Bi,--- ,Bp] stands for a 
block diagonal matrix with matrices Bi, ■ ■ ■ , Bp, and diag[6i, ■ ■ ■ , 6p] is a diagonal matrix with diagonal 
entries 6i, ■ ■ ■ ,6p. The superscripts (•)^, (•)^, (■)*, (■) stand for conjugate transpose, transpose, complex 
conjugate, binary complement, respectively, and V denotes for-all. R+ and C stand for the set of positive 
real numbers and the complex numbers, repectively. dmin is the minimum Euclidean distance between 
two points in the constellation. N and M stand for the number of transmit and receive antennas. 

II. BICMB vv^iTH Constellation Precoding 

Fig. [H represents the structure of BICMB with constellation precoding. First, the code rate Rc = 
kc/uc convolutional encoder, possibly combined with a perforation matrix for a high rate punctured code, 
generates the codeword c from the information bits. Then, the spatial interleaver distributes the coded bits 
into S < min(A^, M) streams, each of which is interleaved by an independent bit-wise interleaver tt. The 
interleaved bits are mapped by Gray encoding onto the symbol sequence X = [xi ■ ■ ■ x^-], where x^. is 
an S* X 1 symbol vector at the /c*^ time instant. In this model, we assume that each stream employs the 
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same modulation scheme. Each entry in the symbol vector belongs to a signal set x C C of size |x| = 2"^, 
such as 2™-QAM, where m is the number of input bits to the Gray encoder. 

The symbol vector x/,. is multiplied by the S x S precoder 0, which is defined as 



where is the P x P unitary constellation precoding matrix that precodes the first P modulated entries 
of the vector x^. When all of the S modulated entries are precoded (P = S), we call the resulting system 
Bit-Interleaved Coded Multiple Beamforming with Full Precoding (BICMB-FP), otherwise, we call it 
Bit-Interleaved Coded Multiple Beamforming with Partial Precoding (BICMB-PP). The symbol generated 
by is multiplied by T which is an 5 x S permutation matrix to define the mapping of the precoded 
and non-precoded symbols onto the predefined subchannels. Let us define bp = [6p(l) ■ • • bp{P)] as a 
vector whose element bp(u) is the subchannel on which the precoded symbols are transmitted, and ordered 
increasingly such that bp{u) < bp{v) for u < v. In the same way, b„ = [&„(!) • ■ • bn{S — P)] is defined 
as an increasingly ordered vector whose element 6„(m) is the subchannel which carries the non-precoded 
symbols. 

The MIMO channel H G C^^^^ is assumed to be quasi-static, Rayleigh, and flat fading, and perfectly 
known to both the transmitter and the receiver. In this channel model, we consider that the channel 
coefficients remain constant for the K symbol duration. The beamforming vectors are determined by the 
SVD of the MIMO channel, i.e., H = UAV^ where U and V are unitary matrices, and A is a diagonal 
matrix whose s*^ diagonal element. As G M^, is a singular value of H in decreasing order. When S 
symbols are transmitted at the same time, then the first S vectors of U and V are chosen to be used as 
beamforming matrices at the receiver and the transmitter, respectively. U and V in Fig. [T] denote the first 
S column vectors of U and V. 

The spatial interleaver arranges the symbol vector x^ as x^ = [x^b^ :x^b„]^ = [^fc,fep(i) ' ' ' ^k,bp(P) '■ 
Xk,b„{i) ■ ■ ■ Xk,b„{s-P)]'^ where x^^bp and Xfc^b„ are the modulated entries to be transmitted on the sub- 
channels specified in bp and b.„, respectively. The 5x1 detected symbol vector = [(r^)"^ : (ffc)^]^ = 
[rk,i ■ ■ ■ rk,p : rk,p+i ■ ■ ■ rk,sY the k*^ time instant is 
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Tk = r0Xfc + rifc 



(2) 



where F is a bloclc diagonal matrix, T = diag[rp, r.„] with diagonal matrices defined as Tp = diag[Afcp(i), 
■ ■ ■ , \(P)], r„ = diag[Afc„(i), ■ ■ ■ , Ab„(5_p)], and rifc = [(n^)^ : (n^)^]^ = [r„,i ■ ■ ■ r„,p : r„,p+i ■ ■ ■ r^^^sf 
is an additive white Gaussian noise vector with zero mean and variance A'^o = N/SNR. H is complex 
Gaussian with zero mean and unit variance, and to make the received signal-to-noise ratio SNR, the total 
transmitted power is scaled as A^. The input-output relation in Q is decomposed into two equations as 

(3) 

The location of the coded bit c^/ within the symbol sequence X is known as k' — > (k,l,i), where k, 
I, and i are the time instant in X, the symbol position in x^, and the bit position on the symbol Xk,i, 
respectively. Let xl denote a subset of x whose labels have 6 G {0, 1} in the i*^ bit position. By using the 
location information and the input-output relation in Q, the receiver calculates the maximum likelihood 
(ML) bit metrics for the coded bit Ck' as 



7''^(rfc,CfcO = min ||rfc - TOxf (4) 
where is a subset of x"^, defined as 

^b' = {x = [xi ■ ■ ■ xsf ■■ Xs\s=i e xl and Xs\s^i G x}- 

In particular, the bit metrics, equivalent to (01) for partial precoding, are represented as 

{mill II - r^exf, if 1 < / < P 
(5) 
min \rk i - A,-a;p, if P + 1 < I < S 

where V^^'* is a set which is mapped from the set by a surjective function /(x), for x = [xi ■ ■ ■ x^]^, 
defined as 

/ (x) = [xi ■■■ Xpf, 

and / is an entry in b„, corresponding to the subchannel mapped by T. Finally, the ML decoder makes 
decisions according to the rule 



argmm V'7''*(rfc,Cfc/). (6) 
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III. Diversity Analysis 

Since BER in BICMB is bounded by the union of the PEP corresponding to each error event flU, the 
calculation of each PEP is needed. In particular, the overall diversity order is dominated by the pairwise 
errors which have the smallest exponent of signal-to-noise ratio in PEP representation. In this section, we 
calculate the upper bound to each PEP corresponding to the pairwise errors. 

A. BICMB with Full Precoding 

Based on the bit metrics in (HI), the instantaneous PEP between the transmitted codeword c and the 
decoded codeword c is calculated as 



Pr(c c|H) = Pr ( ^ min ||rfe - rGxf > ^ 

\ k' ^^^i' k> 



min \\Yk - rOxf (7) 



where Ck' and ca,/ is the coded bit of c and c, respectively. We define du as the Hamming distance between 
c and c. It is assumed that the dn coded bits are interleaved such that they are placed in distinct symbols. 
In addition, we know that the bit metrics corresponding to the same coded bits between the pairwise 
errors are the same. Based on the assumption and the knowledge, © is re-written as 



Pr(c c|H) = Pr ( V min ||rfc - r0xf > V min ||rfc - r0xf 

xk^H ^^^"k' k,d„ ^^^Sfc/ y 



(8) 



where Xlfe djj stands for the summation of the dn values that correspond to the different coded bits between 
the codewords. 



Let us define x^ and x^ as 



Xfc = arg min ||rA; — r0x||^ 

(9) 

Xfc = arg min \\rk — r0x||^ 



where Ck' is the complement of Ck' in binary codes. It is easily found that x^ is different from x^ since 
the sets that the symbols belong to are disjoint, as can be seen from the definition of In the same 
manner, we see that x^. is different from x^. With x^. and x^, we get the following expression from ([8]) as 



Pr (c ^ c|H) = Pr I J] ||r,. - r0i,f > J] ||r,. - r0x,fj . (10) 



Based on the fact that Hr^ — rOx^H^ > ||rfc — r0Xytp and the relation in equation (flOl ) is upper- 
bounded by 



Pr(c ^ c|H) < Pr ( /3 > ^ ||r0(xfc 

\ k,dH 



(11) 



where 



/3 = - ^(xfc - Xfc)^e^rn, + nf r0(x, - X,). 



k,dH 



Since /3 is a zero mean Gaussian random variable with variance '^NoJ^kdn 11-'^® (^fc ~ CH) 1^ 

replaced by the Q function as 



Pr(c c|H) < Q 



\ 



E ||re(x»-xt)||2 

k,d[j 

2iVo 



v 
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The numerator in (|T2|) is rewritten as 

s s 

lire(x, - X,) f = E xl\d,,s\' = E E 



(13) 



kAi 



kAji s=l 



s=l kAi 



where = [d^^i ■ ■ ■ d^^sV — ®(xa: — x^). Using an upper bound to the Q function, we calculate the 



average PEP as 



Prfc ^c) < E 



exp 



s=l k,du 



4iVn 



V 



/J 



(14) 



In im, we have shown that equations with such form as (fT4l) have a closed form expression of an upper 
bound. We provide a formal description below. 

Theorem 1: Consider the largest S < min(A^, M) eigenvalues jjs of the uncorrected central M x N 
Wishart matrix that are sorted in decreasing order, and a weight vector = [0i ■ ■ ■ 0^]^ with non- 
negative real elements. In the high signal-to-noise ratio regime, an upper bound for the expression 
£^[exp(— 7 Ef=i 4>sf^s)] which is used in the diversity analysis of a number of MIMO systems is 

s 



E 



exp -7 0s/is 



s=l 



where 7 is signal-to-noise ratio, C is a constant. 



riin = min{0i, ■ ■ ■ , (f)s}, and 5 is the index to the first 



non-zero element in the weight vector. 

Proof: See m. □ 
By calculating the weight vector whose s*^ element is Ylikdu M^.sl^' evaluate the diversity order of a 
given system. In particular, if the constellation precoder is designed such that 

|4,i|' = |6>r(xfc - Xk)P > 0, V(xfe,x,.) (15) 

where d\ is the first row vector of the unitary precoding matrix 0, we see that \dk,i\^ > 0, resulting 

in the full diversity order of NM. Therefore, (fT5l) is a sufficient condition for the full diversity order of 
BICMB-FP. 

B. BICMB with Partial Precoding 

The bit metrics in ^ lead to the PEP calculation as 

Pr(c ^ c|H) = Pr(ri > Ta) (16) 

where 

Ti = min ||r^ — rp0x||^ + min |rfc^; — A^-x|^ 

Kd„^^ ''^'^4' k,dH,n ^^^"'fe' 



T2= ^ min ||r^ - rp0x||^ + ^ min \rk,i-\ix\'^ 

dn ' ^fc dH st^'^d for the summation over the dH,p and d^,™ bit metrics corresponding to the 
different coded bits carried on the subchannels in bp and b„, respectively. By using the appropriate system 
input-output relations, the PEP is written as 

Pr(c ^ c|H) = Pr(/3 > k) (17) 

where /5 = /3p + 

(^P = - Yl '^^'^'bp ~ Xfe,bj^0^rpn^ + (n^)^ rp0(xfe,bp - Xfc,bJ, 

k,dH,p 
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Pn = - >^iixk,i - Xk,i)*nk,i + \i{xk,i - Xk,i)nl^i, 

k,dH,n 

and 

K = ^ \\Tp@ (xfc,bp - Xfc,bJ f + ^ |A[ (xfc,/ - Xk,i) 1^. 

Since P in (fTTl ) is a Gaussian random variable with zero mean and variance 2Nqk, the PEP can be 
expressed in a similar way as (fT2l) with the Q-function. In addition, if we define a as 

p S-P 

= '^ftpM Mfc,rP + C?Ln '^b„{r)(^b„ir) (18) 

where = [dk,i ■ ■ ■ dk,pV = ® (^fc.bp — x^ bp), and is the number of times the s^^ subchannel is 
used corresponding to „ bits under considertion, then we can see that o > k. Finally, the average PEP 
is calculated as 



Vxic^c) < E 



a 



2 V 4M 



(19) 



To determine the diversity order from a, we need to find the index to indicate the first non-zero element 
in an ordered composite vector which consists of J^kdn \^k,r\'^ and ab„(^r) as in Theorem [TJ If dH,p = 0, 
the first summation part of a vanishes. In this case, the first index is 

S = mm{s : a, > for s G {6n(l), ■ ■ ■ , &„(^ - P)}}. (20) 

In the other case of rf^/p > 0, we see that bp and X/t^bp are obviously different for the same reason 
as in the previous section. If the constellation precoder satisfies the sufficient condition of (fT5l) . the term 
with always exists in a. Therefore, 5 for the case of dn^p > is 5 = min(6p(l), 5') where 5' is 

obtained in the same way as (|20l ). 

Example of Determining Diversity Order : 

In this example, we employ 4-state 1/2-rate convolutional code with generator polynomials (5,7) in 
octal representation mN = M = S = 3 system. Two types of spatial interleavers are used to demonstrate 
the different results of the diversity order. A generalized transfer function of BICMB with the specific 
spatial interleaver and convolutional code provides the a-vectors for all of the pairwise errors, whose 
element indicates the number of times the stream is used for the erroneous bits [8J. In particular, due to 
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the fact that dH,p = J2r=i '^bp(r) and dH,n = J2r=i '^bn{r) where is the s*^ element of the a-vector, the 
generalized transfer function is also useful in the analysis of BICMB-PP. Hence, we rewrite the transfer 
functions of the systems from [[8]|, where a, b, and c are the symbolic representation of the I''*, 2"°', 3'"'^ 
stream. The spatial interleaver used in Ti is a simple rotating switch on 3 streams. For T2, the u*^ coded 
bit is interleaved into the stream Sniod(«-i,i8)+i where si = ■ • • = sq = 1, st = ■ ■ ■ = S12 = 2, S13 = ■ ■ ■ = 
S18 = 3 and mod is the modulo operation. Each term represents the «-vector, and the powers of a, b, c 
indicate the elements of a-vector. 

Ti = Z'{a%^c + a%c^ + ab^c^) + Z\a%''c + ■ ■ ■ ) 

+ Z\2aH^c + 2aH'^c^ + 2a%^c^ + ■ ■ ■ ) (21) 
+ Z^{a^b^ + b\^ + a^c^ + ■■■) + ■•• 



T2 = Z^(a^ + aV + 0^6^+ 

b^ + a^c^ + b^c^ + a^e" + b'^e' + c^) 
+ Z\a%'' + ?>a%^ + + a^c^ + ?>a%''c^+ (22) 
b^c^ + ?,a^c^ + ?>b^^ + + b\^) + ■ ■ ■ 

Consider the case bp = [12]. We see that all of a-vectors of Ti show dn^p > 0, leading to 5 = 1. 
Therefore, the diversity order of the Ti BICMB-PP system with bp = [1 2] achieves the full diversity order 
while BICMB without constellation precoding |l8l, or PPMB without bit-interleaved coded modulation 
(BICM) loses the full diversity order |l9l [fTOll . However, T2 has [0 5] which shows dH,p = 0, resulting 
in S = 3. Therefore, the diversity order of the T2 BICMB-PP system with bp = [1 2] does not achieve 
the full diversity order. 

The same analysis for bp = [1 3] results in the diversity order of 9, and [2 3] results in 4 for the transfer 
function Ti. Similarly, both of [13] and [2 3] result in the diversity of 4 for T2. As a consequence, we 
find that proper selection of the subchannels for precoding, as well as the appropriate pattern of the spatial 
interleaver, is important to achieve the full diversity order of BICMB-PP. 

IV. Simulation Results 

Monte-Carlo simulations were performed to verify the diversity analysis in Section Hill Throughout the 
simulations, we used the precoding matrices in |l9l, [fTOll which meet the sufficient condition to achieve 
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the full diversity order of ([T?] ). Fig. [2l depicts the simulation result for 2 x 2, 3 x 3, and 4x4 BICMB 
and BICMB-FP with 64-state convolutional code punctured from 1/2-rate mother code with generator 
polynomials (133, 171) in octal representation. In [8], we showed the maximum achievable diversity 
order of BICMB with an i^^-rate convolutional code is {N - \S ■ Rc] + 1)(M - \S ■ R^] +1). In this 
example, the maximum achievable diversity order of the three BICMB systems is 1. However, Fig. |2] 
shows that BICMB-FP achieves the full diversity order for any code rate. Fig. [3] depicts the simulation 
results of BICMB-PP given in the example of Section IIII-B[ The diversity orders of the BICMB systems, 
Ti and T2 are 4 and 1, respectively. We see that the simulation results match the analysis in IIII-B[ 

To compare coding gain between BICMB and BICMB-PP that achieve the full diversity order, we show 
in Fig. m the BER performances with the 4x45 = 4 systems. The used generator polynomials of 64- 
state, 1/2 and 1/4-rate convolutional codes are (133, 171) and (117, 127, 155, 171) in octal representation, 
respectively [fTTll . As shown in the figure, BICMB combined with constellation precoding shows larger 
coding gain with a large number of antennas and at the higher transmission rate. 

V. Conclusion 

We investigated the diversity order of BICMB combined with the constellation precoding scheme, by 
calculating pairwise error probability. Using the analysis, we presented the resulting diversity order of 
the given examples. The analysis can be used to determine the precoding configuration from the given 
BICMB implementation to get the full diversity order. We provided simulation results that proves the 
analysis. In addition, the simulation showed that BICMB-PP outperforms BICMB with a large number 
of antennas and at the higher transmission rate. 
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Fig. 1. Structure of Bit-Interleaved Coded Multiple Beamforming with Constellation Preceding. 
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Fig. 2. BER comparison between BICMB and BICMB-FP with 16-QAM, and 64-state punctured convolutional code. 
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Fig. 3. BER vs. SNR for BICMB-PP with 3x35 = 3, 4-QAM, and 4-state 1/2-rate convolutional code. 
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Fig. 4. BER comparison between 4x45 = 4 BICMB and BICMB-PP at 4 and 8 bits/channel use. 



